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ABSTRACT 

We show that the proper motion of the Becklin-Neugebauer (BN) object is consistent with its dy- 
namical ejection from the 9 1 Ori C binary, contrary to recent claims by Gomez et al. Continued radio 
observations of BN and future precise astrometric observations of 9 1 Ori C with SIM and the Orion 
Nebula Cluster with GAIA can constrain the properties of this ejection event, with implications for 
theories of how the nearest example of massive star formation is proceeding. 

Subject headings: stars: formation — stars: kinematics 
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1. INTRODUCTION 

Understanding massive star formation remains one of the most challenging and important problems of contemporary 

astrophysics (Beuther et al. 2007; Zinnecker & Yorke 2007). The complexity of the process means that massive star 

formation theories, such as the turbulent core model (McKee & Tan 2003), the competitive accretion model (Bonnell 

& Bate 2006) and stellar coalescence model (Bonnell et al. 1998; Clarke & Bonnell 2008) require close testing against 

observed systems. The closest forming (i.e. accreting) massive star is thought to be radio source I (Menten & Reid 1995) 

O ■ within the Orion Nebula Cluster (ONC), at a distance of 414 ± 7 pc (Menten et al. 2007, adopted throughout), in the 

,H | Kleinmann-Low (KL) region. As reviewed by Tan (2008), this source has been used as observational evidence in support of 

all three of the above theories. Part of this confusion is due to the Becklin-Neugebauer (BN) object, 9.9"to the NW (Fig. 1), 

i ^, ' which is a fast moving (radio- ONC-frame proper motion of /ibn = 13.2 ± 1.1 masyr -1 , i.e. W2D.BN = 25.9 ± 2.2 kms -1 

I towards P.A. BN = -27°. 5 ± 4°, Plambeck et al. 1995; Gomez et al. 2008) embedded B star (L B n = (2.1 - 8.5) x 10 3 i Q , 

£NJ ■ Gezari, Backman & Werner 1998, equivalent to a zero age main sequence mass TOBN,zams = 9.3 ± 2. CM©). This proper 

^ ' motion implies that BN has been moving through the KL region and made a close, possibly coincident, passage with 

~* source I about 500 years ago. Thus to understand the nearest example of massive star formation, we need to understand 

' ■ the origin of BN's motion. 

Including the (+21) — (+8) = +13 km s _1 radial velocity of BN with respect to the ONC mean (Scoville et al. 1983; 
**! ! Walker 1983), the 3D ONC-frame velocity of BN is u 3D ,bn = 29 + 3 km s _1 , and its kinetic energy is E BN = (8.3 ± 2.3) x 
!"***■ ' 10 46 (mBN/10M Q ) ergs. BN is very likely to have formed somewhere in the ONC and then attained its high speed by a 
CD close interaction with a massive multiple stellar system followed by dynamical ejection (Poveda, Ruiz & Allen 1967). 

Tan (2004) proposed BN was launched from the 9 1 Ori C binary (also shown in Fig. [TJ, since this is the only stellar 
system in the ONC known to have all the physical properties required by this scenario: (1) a location along BN's past 
trajectory (§5); (2) an (optical)-ONC-frame proper motion (fJ,e 1 c = 2.3 ± 0.2 masyr" 1 , van Altena et al. 1988, i.e. 
i^d^c — 4-5 ± 0.4 km s _1 , towards P.A.giQ = 142°. 4 ± 4°) that is in the opposite direction to BN (the direction to BN 
from 9 1 Ori C is a P.A.= —30°. 949) and is of the appropriate magnitude (the dynamical mass of BN implied by this motion 
jrt ■ agrees with the estimate of mBN.zams and is maN.dyn = 8.6 ± I.OMq assuming negligible error in m e i c — 49. 5M© and 
negligible motion of the pre-ejection triple system in this direction; a pre-ejection motion of 0.35 mas/yr along this axis (321 
would contribute an additional 1.5M Q uncertainty); (3) primary (mgic-i = 34Mq) and secondary (mgi c _ 2 = 15.5Mq) 
masses greater than tobn (Kraus et al. 2007); (4) a semi-major axis of a — 17.0 + 5.8 AU (Patience et al. 2008) and thus a 
total orbital energy (E tot = Gmeic*_im 9 i c _ 2 /(2a) = (2.7 + 0.9) x 10 47 ergs) greater than the sum of BN's kinetic energy 
and 9 1 Ori C's kinetic energy (1.00 x 10 46 ergs) (see Tan 2008 for a review). Note, 9 1 Ori C's recoil in this scenario is 
large enough to remove it from the Trapezium region (see Pflamm-Altenburg & Kroupa 2006 for theoretical studies of the 
dynamical decay of Trapezium-like systems) and may be enough to eject it from the ONC completely, with implications 
for the effectiveness of its ionizing feedback on disrupting the star cluster formation process. 

Rodriguez et al. (2005) and Bally & Zinnecker (2005) proposed BN was launched from an interaction with radio source 
I, which would require this system to be a massive binary, recoiling away from any large scale (> 100 AU) gas that it was 
originally accreting. Gomez et al. (2008) used the relative motion to BN with respect to source I to claim that BN could 
not have made a close passage with 9 1 Ori C, excluding this possibility at the 5-10 a level. 

We show in [J5]that if BN's motion is considered in the reference frame of the ONC, then a close (coincident) passage 
with 9 1 Ori C is allowed by the data, which permits the scenario of dynamical ejection of BN from 9 1 Ori C. In ffJIwe 
discuss the potential of future high precision astrometric measurements to constrain the properties of BN's dynamical 
ejection, which then constrain BN's interaction distance with source 7, the mass of source I, and thus the strength of tidal 
perturbations on the massive protostar during this encounter. 
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2. ASTROMETRY OF BN IN THE ORION NEBULA CLUSTER 
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Fig. 1. — This diagram shows the positions of the Trapezium stars 9 1 Ori A, 6 1 Ori B, 9 1 Ori C and 8 1 Ori D that make up the core of the 
ONC. The positions of radio sources I and BN are also shown. The coordinates are relative to the present position of source I (a(J2000)=05 
35 14.5141, 5(J2000)=-05 22 30.556) (Gomez et al. 2008). The proper motions relative to the cluster of BN (Gomez et al. 2008) and 9 1 Ori C 
(van Altena et al. 1988) are indicated with the arrows. Past trajectories (dashed line) and la uncertainties (dotted lines) are drawn. 



To determine BN's past trajectory through the ONC we use the absolute proper motion of BN (fi a cosS — —5.3 ± 
0.9masyr _1 , fig — 9.4 ± l.lmasyr -1 (la errors); Gomez et al. 2008) and then correct for the motion of the ONC (mean 
of 35 radio sources within central 0.1 pc of ONC: fi a cosS — +0.8 ± 0.2masyr _1 , fig — —2.3 ± 0.2masyr _1 ; Gomez et al. 
2005). The ONC-frame proper motions are shown in Fig. [TJ One sees that the past trajectory of BN through the ONC 
overlaps within the la errors with the present position of 8 1 Ori C. Given the motions of BN and 9 1 Ori C, the time of 
coincidence (i.e. when the dynamical ejection took place) was 4530 years ago, i.e. about 174 orbital periods of 8 1 Ori C 
(although the orbital period is only poorly constrained at present to 26 ± 13 years, Patience et al. 2008). 

Gomez et al. (2008) excluded a coincidence between BN and 9 1 Ori C because they used the motion of BN with 
respect to source / (which is measured using relative astrometry to greater accuracy so has smaller error bars), but did 
not allow for the fact that their data indicate that source / is moving. In the ONC frame this motion is claimed to be 
H a cosd = — 3.7 ± 1.2masyr _1 , /j,g = —3.4 ± 1.3masyr _1 , corresponding to /ir = 5.0 ± 1.3 masyr" 1 (i.e. 9.9 ± 2.6km s _1 ) 
towards a P. A.= +133° ± 16°. 

We note, as a separate point, that source lis elongated along the NW-SE axis, i.e. towards P.A.~ +135° (Reid et al. 
2007). If the source exhibits variability affecting the location of the centroid of its emission, then this could lead to an 
apparent, but false, proper motion. This effect is a potential source of additional uncertainty in the motion reported for 
source / (and for source n) by Gomez et al. (2008). 

Source / is thought to be a massive protostar and a large proper motion would be interesting for theories of massive 
star formation. Furesz et al. (2008) measured the distribution of radial velocities in the ONC, finding it could be well 
fit by a Gaussian with auj = 3.1 kms -1 , for both the entire cluster and for stars within a 15' radius of the Trapezium. 
Assuming an isotropic velocity distribution, the proper motions should exhibit a Gaussian distribution of motions with 
o"2D = 4.4 km s _ . In comparison, Source Fs claimed motion of 9.9 ± 2.6km s _1 is (2.3 ± 0. 6)021), i-e. not significantly 
larger than expected of a typical cluster member. Note, Jones & Walker (1988) found a 2 D — 2.9 km s _1 from direct 



observation of proper motions (adjusted to rfoNC = 414 pc), for which source Ps motion would then be (3.4 ± 0.9)(T2D. 
Gomez et al. (2005) found <72D — 7.6 km s~ based on proper motions of 35 radio sources, for which source Fs motion 
would then be (1.3 ± 0.3)(T2d- We conclude, in contrast to Gomez et al. (2008), that it is premature to claim that source 
I has an anomalously large motion compared to other ONC stars. 

3. POTENTIAL OF HIGH PRECISION ASTROMETRY WITH SIM 

For wide angle absolute astrometry, SIM should be able to achieve a parallax accuracy of about 5 /j,a.s. Assuming a 
distance of about 400 pc, this will allow a parallax distance measurement accurate to 0.2%, i.e. 0.9 pc. 

Once the motions of the primary and secondary components of 6 l Ori C due to their binary orbit are accounted for, 
then the absolute proper motion of the system should be known to an accuracy of a few /ias/yr. By averaging over 
many stars, an even greater accuracy should be achievable for the absolute proper motion of the ONC with GAIA. Since 
9 1 Ori C is moving at a few mas/yr in the ONC frame (van Altena et al. 1988), then the accuracy of the position angle 
of the direction of motion would be ~ 0.06°. Presently it is only known to about 4°. 

If, as seems very likely, BN was ejected from 9 l Ori C, it should have been ejected in exactly the opposite direction 
to 9 1 Ori C's motion as measured in the center of mass frame of the pre-ejection triple system. Comparison of the 
ONC-frame motion of 9 l Ori C with the present position and ONC-frame motion of BN, will yield information on motion 
of the pre-ejection triple system and any accelerations experienced by the stars since ejection. 

The expected size of pre-ejection triple system proper motion is uncertain. If the system (with total mass ~ 60M Q ) was 
in kinetic energy equilibrium with the other ONC stars (with, say, typical mass I.OMq and 02D = 4.0km s~ ), then we 
would expect it to have a plane of sky motion ~ 0.52 km s~ equivalent to a proper motion of 0.26 mas/yr. The observed 
proper motion dispersion of bright (V < 12.5), i.e. massive, stars is 0.70±0.06 mas/yr (van Altena et al. 1988). Assuming 
a 0.5 mas/yr proper motion for the pre-ejection triple system, of which 0.35 mas/yr would be expected to be tangential 
to the ejection axis, implies that the ONC-frame proper motion vectors of 9 1 Ori C and BN would be misaligned by 10° 
from direct opposition. The current observed misalignment is 10° ± 6°. Thus, in the limit that subsequent accelerations 
are negligible, high precision ONC-frame proper motions of 9 1 Ori C and BN (the latter expected from continued radio 
observations) can constrain the motion of the pre-ejection triple system. 

The expected gravitational accelerations of 9 1 Ori C and BN depend on the distribution of mass in their surroundings. 
Their trajectories are taking them away from the ONC center, so they will be experiencing a deceleration associated with 
climbing out of the cluster potential. This effect is largest for BN, but it is still small. BN has moved 0.12 pc (projected) 
from the ejection site, and if the enclosed mass is 500 M Q (likely to be a conservative upper limit, e.g. Hillenbrand & 
Hartmann 1998), then for a starting velocity of 30 km s _1 , it would have decelerated by only 0.6 km s _1 . 

Close passage with individual stars can also cause more significant accelerations. 9 1 Ori C's trajectory may have brought 
it into relatively close proximity with 9 1 Ori A (a B0 star, i.e. 16M Q , 13" to the NW, with a visual companion at 100 AU 
of 4Mq and a spectroscopic companion at ~ 1 AU of ~ 3M Q , Schertl et al. 2003). However, the relative motion of these 
stars is only about 1.2 mas yr -1 (van Altena et al. 1988) so that the time of closest approach would have been about 
10 4 yr ago, long before the proposed interaction of 9 1 Ori C with BN. 

More importantly, BN made a close passage to source I about 500 years ago. From the bolomctric luminosity of 
the KL region, source I is estimated to have a protostellar mass of about 20 Mq. As an example of the magnitude 
of the deflections that can be expected, treating BN as a massless test particle, its deflection angle due to source I is 
2.25°(m/ j */2OAf0)(6/lOOOAU)~ 1 (vBN/3Okms~ 1 )~ 2 , where b is the initial impact parameter and Vbn is the velocity of BN 
relative to source I. A direct trajectory from 9 1 Ori C's present position (ideally this would be measured from 9 1 Ori C's 
position at the time of ejection) to BN's position has a closest projected separation from source Fs present position of 
1.5" (about 600 AU). Thus an accurate astrometric solution of this system presents us with the unique opportunity of 
constraining the dynamical mass of source 7, the nearest massive protostar, in combination with the true (unprojectcd) 
distance of closest approach. The true distance of closest approach is important for evaluating the tidal effects of BN on 
source Ps accretion disk, which are likely to have enhanced accretion to the star (Ostriker 1995; Moeckel & Bally 2006). 
Such enhanced accretion is likely to have led to enhanced protostellar outflow activity, thus explaining the ~ 1000 yr 
timescale of the "explosive" outflow from this region (Allen & Burton 1993; Tan 2004). 

4. CONCLUSIONS 

We have reviewed the latest evidence that BN was dynamically ejected from the 9 1 Ori C binary, finding that 9 1 Ori C 
has all the physical properties expected in this scenario. We showed that the trajectory of BN is also consistent with this 
scenario, in contrast to recent claims by Gomez et al. (2008). We discussed how high precision astrometry of 9 1 Ori C 
with SIM can yield information on the pre-ejection velocity of the system and the size of any subsequent deflections, in 
particular that of BN caused by close passage with source /, the nearest massive protostar. 

JCT acknowledges support from NSF CAREER grant AST-0645412 and a grant from NASA for SIM Science Studies. 
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